Scaling of entanglement entropy and superselection rules 
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Particle number conservation in fermionic systems restricts the allowed local operations on bi- 
partite systems. We show how this restriction is related to measurement entropy of particle fluctua- 
tions and compute it for several regimes of practical relevance. The accessible entanglement entropy 
restricted by particle number conservation is equal, to leading order, to the full entanglement en- 
tropy. The correction is bounded by the log of the variance of particle number fluctuations. The 
results are applied to generic Fermi systems in dimension d > 1 as well as to several critical systems 
in d = 1. 



Entanglement entropy is a quantity which is often used 
in quantum information theory to characterize the de- 
gree of entanglement between different parts of a quan- 
tum system that share the same microscopic variables [l[ . 
Prior to its adoption in the field of quantum information, 
entanglement entropy has been introduced and studied in 
field theory as a possible contribution to the entropy of 
black holes @, H, H| , where it was called "geometric en- 
tropy." Recently it emerged as a quantity of interest in 
many-body theory, where it was employed as a measure 
of large-scale, nonlocal correlations signaling criticality 
(for a review of these developments see Ref.0] and 
references therein). 

The non-local correlations between entangled quantum 
variables are revealed by measurements made locally on 
two (or more) subsystems of a quantum system. The 
locality of observables used to detect entanglement is 
of particular importance in applications such as telepor- 
tation [3] and quantum cryptography Q, which require 
that the two subsystems do not interact after the initial 
state is prepared. However, our ability to probe entan- 
glement by a combination of local measurements is often 
constrained by conservation laws (e.g. of particle num- 
ber, charge or spin), rendering some local observables 
physically inaccessible. 

Constraints on the local operations (known as super- 
selection rules) can limit experimentally accessible entan- 
glement. The interplay of super-selection and entangle- 
ment was first analyzed in Ref.0], where the accessible 
entanglement was quantified by averaging of entangle- 
ment entropy over super-selection sectors. The quantity 
defined in Ref.@l, which we refer to as accessible entan- 



glement entropy was further discussed in Refs. 
A different point of view on accessible entanglement was 
put forward in Ref. [l]| , where super-selection constraints 
on local operations are treated as a resource that can be 
employed for hiding information in correlations which are 
blocked from local probing. 

To understand the general relation between accessi- 
ble entanglement entropy and the full entanglement en- 
tropy, in this Letter we present general results for fermion 
systems, interacting or noninteracting. Fermion systems 



play a special role in the theory of entanglement entropy 
because, on one hand, they are many-body systems ex- 
hibiting rich and interesting behavior and, on the other 
hand, they provide a good model of experimentally rele- 
vant settings. In particular, the scaling of the full entan- 
glement entropy in Fermi systems depends on the nature 
of the state. If the system is gapped, the entropy scales 
with the area of the boundary of the region, S oc L d ~ x , 
where L is the linear size of the region (c.f. [3, [H|)- 
However if the system is in a gapless (metallic) state, the 
entropy scales as S cx L d ~ 1 logL [13] • Furthermore, 
dynamics of fermionic systems can be used to generate 
complex entangled states. In particular, generation of en- 
tanglement in mesoscopic conductors occurs naturally as 
a result of elastic scattering of electrons on barriers and 
disorder potential 
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19j . A scheme for detecting elec- 
tronic entanglement which can discriminate between oc- 
cupation number entanglement and mode entanglement 
has been proposed in [201 ]. 

One particularly attractive aspect of Fermi systems is 
that the full many body entanglement entropy of free 
fcrmions can be linked to experimentally accessible quan- 
tities, such as particle number fluctuations [2l[ or current 
fluctuations |22j . The latter quantity was discussed in de- 
tail in a recent proposal [22I ] of an experiment to measure 
the entanglement entropy by detecting electric noise gen- 
erated in a process of connecting two Fermi seas through 
a quantum point contact (QPC). 

While proposals such as [U, [52] as well as studies of 
scaling of entropy in fermion systems 3, 13] are of in- 
terest from the quantum information perspective, they 
do not make a distinction between entanglement entropy 
and the accessible entanglement entropy. Given that only 
the accessible entropy can be used as a source of entangle- 
ment for quantum information applications, here we set 
out to investigate the effect of super-selection rules on 
the scaling of the accessible entanglement in many-body 
systems of fermions. We find that, under very general 
assumptions, super-selection rules result in a difference 
between the accessible and full entanglement entropies 
which is small in a relative sense, provided that the num- 
ber of Fermi particles contributing to the entropy is large, 
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N ;§> 1. It means that, in essence, the schemes pro- 
posed for measuring the full entanglement entropy can 
also yield a good estimate of the accessible entanglement 
entropy, and vice versa. 

To state our results in a quantitative form, let us recall 
that entanglement entropy is defined for a system parti- 
tioned into two parts A and B. The quantum state of the 
system is projected on A, giving the reduced density ma- 
trix pa — Trap, where all degrees of freedom outside A 
have been integrated out. Entanglement entropy is then 
given by the von Neumann formula: Sa = — Trp^ log pa- 
When super-selection rules are present, a natural quan- 
tity to consider is Sa s [S|, which is obtained by averag- 
ing entanglement entropy over subspaces with fixed con- 
served quantity (in our case, particle number) as will be 
described below. We show that if the state p has a fixed 
number of particles then 



S A - AS < ST < S A , 



A5=ilog [2ne(C 2 + ±)} (1) 



where C2 = ((Na— (Na)) 2 } is the variance in the number 
of particles in subsystem A. An immediate consequence 
is that since typically C2 can not grow faster than polyno- 
mially with the size of the subregion, the correction AS 
to the entropy due to particle conservation is at most 
logarithmic in the volume. 

In many-body systems, where entropy scaling is of in- 
terest, entanglement entropy Sa of translationally invari- 
ant systems (gapless as well as gapped) in space dimen- 
sion d > 1 typically scales at least as the boundary area 
where L is the size of the region A. Given that 
AS in Eq.([T]) grows at most as logL, we see that the 
correction to the entropy from the terms violating super- 
selection is sub-leading to the entropy Sa- The situation 
may be more complicated in dimension d = 1, where for 
critical systems, described by a conformal field theory, 
entropy of a region of length L scales as log L [3(| . Still, 
even in this case, if the fermion density-density correla- 
tions decay as 1/r 2 or faster, the quantity logC2 grows 
at most as log log L, and thus AS* is again sub-leading to 
Sa- Several systems of interest exhibiting this behavior 
are analyzed below. 

Perhaps surprisingly, we found that S r ^ s may be com- 
puted explicitly in many experimentally relevant settings. 
This is largely due to the fact that the difference Sa—Sa 8 
is nothing but the measurement entropy S m of particle 
number in subregion A. We analyze the accessible en- 
tanglement entropy for several cases: Entanglement gen- 
erated when two Fermi seas are connected via a QPC 
for a time At and then disconnected, entanglement gen- 
erated by a dc current in a QPC biased by voltage V 
during time At, entanglement in a Luttinger liquid and 
in a d-dimensional free fermion system. The results for 
S m = Sa — S™ s are summarized in Table I. 



System 


C C nres 
Om — OA — iJA 


Fermi seas connected via a QPC (Fig. la) 
DC voltage-biased QPC (Fig. lb) 
Chiral Luttinger liquid 
Fermions in dimension d 


Ioglog(At/r) 
log(AteV/fc) 
log log(Z) 
log(L d_1 log L) 



TABLE I: The difference between full and accessible entropy 
for the systems described in the text (r is a short time cutoff, 
L — kpL, where Jcf is Fermi momentum and L is system 
size). 
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FIG. 1: Schematic many-body evolution in a quantum point 
contact (QPC), generating entanglement and current fluctu- 
ations. Two fermionic reservoirs are coupled through a tun- 
able tunnel barrier, switching between fully closed and fully 
open states (a), or being time- independent (b). In the case 
(a), the left and right reservoirs, which are connected during 
< t < At, and disconnected at earlier and later times, are 
maintained at equal chemical potentials. In the case (b) the 
QPC is biased by a DC voltage V, with transmission fixed at 
a constant value < D < 1. 



SUPERSELECTION AND ENTANGLEMENT 
FOR TWO FERMIONS 



We begin with some general comments on the nature of 
entanglement of fermions and the effect of particle num- 
ber conservation. To motivate the definition we use for 
the accessible entanglement entropy, we briefly discuss 
what entanglement of fermions means, and how conser- 
vation laws enter the picture. In that, for reader's sake, 
we restate some of the observations made in Ref. [s| . 

Fermionic wavefunctions are always antisymmetrized, 
but this property by itself does not imply entangle- 
ment. To illustrate the role of antisymmctrization and 
of particle conservation we focus on the simplest case of 
just two fermions. Consider a one dimensional interval 
< x < L, partitioned into equal halves A and B. We 
define A(x) = lifO<x<L/2 and A = otherwise and 
similarly B(x) = 1 ifL/2<a;<_Lan.B = otherwise. 
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Let us compare the properties of two antisymmetric 
wavefunctions, 



^(x,y) = ^(A(x)B(y)-A(y)B(x)) 



and 



V2L 



(e 1 



(2) 



(3) 



The state tpi is obtained by putting one fermion in box A 
and another in box B. While this state may appear "en- 
tangled" due to the antisymmetrization, it is of course 
not entangled. The reduced density matrix in B is just 
a projection on a single particle with wavcfunction B(x), 
and thus it is a pure state. In other words, lack of in- 
formation about the state in A, does not increase the 
entropy of the state in B. Thus, anti-symmetrization of 
the wavefunction by itself does not produce entanglement. 

In contrast, the state ip2, obtained by putting one par- 
ticle in the state -h= and another particle in the state 



e l 



contains entanglement. To see this we decom- 



pose ip2 as follows: 



ip2(x, y) = F AA + F B b + F AB + F B a 



(4) 



where 



Faa = vfe (e ^ -e^)A(x)A(y) (5) 
F B B = ^ z (e^-e^)B(x)B(y) 
Fab = ^(e^A(x)B(y) - e^A(y)B(x)) 
F BA = ^ Z (^B(x)A(y) - e^B(y)A(x)) 

The Faa and Fbb wavefunctions represent the possi- 
bility that both fermions are simultaneously in A or in B. 
Let us now measure the number of particles in A. If the 
result of this measurement is or 2, the wavefunction is 
collapsed into F A a or Fbb- However, if we find just one 
particle in A, we have collapsed the wavefunction into 
the state Fab + Fba ■ 

The wavefunction Fab + Fba is maximally entangled, 
in complete analogy with a singlet spin state of Einstein- 
Podolsky- Rosen form. Indeed, we can denote | X)A and 

| J,)^ as having a particle in A in the modes y^e 1 ^ A(.t) 

and iJ~^A(x), and similarly for | ])b and | [)b (the form 

e^r^ was chosen so that a{] \ [)a = 0). In this notation 
we have 

F A B+F B A = ^{\])A\i)B + \i)A\])B) (6) 

Locally measuring whether the "up" or "down" mode is 
occupied, and locally moving a particle between these 
states (for example by coupling to a potential) are valid 



physical operations corresponding in spin notation to o~ z 
and o~ x . 

This entanglement can be contrasted with the proper- 
ties of the two-particle state Faa + Fbb- While the state 
Faa + Fbb is formally entangled, this entanglement can 
not be revealed by Bell measurements. The states anal- 
ogous to I and f here will be "having no particles in 
A" and "having two particles in A" , which could be also 
associated with spin variables o~ z = ±1 and directly mea- 
sured. However, because of particle number conserva- 
tion, we are blocked from rotating between these "up" 
and "down" states, and so u x is inaccessible. This blocks 
us from a necessary ingredient to perform Bell measure- 
ment, rendering this entanglement "unuseful" . We have 
thus seen that in this example, with probability 1/2, we 
can perform measurements that violate Bell inequalities 
on %p 2 - 

Consider a large number n of systems with the wave- 
function ip2- Let p denote the probability of measuring 
one charge in A in state tp2 ■ If we now measure the num- 
ber of particles in A, we will average pn cases of (maxi- 
mally entangled pairs) F A b + Fba- This motivates defin- 
ing the entanglement of the system by averaging over the 
probabilities to fall in different super-selection sectors. 

The discussion above, was concerned with the case of 
copies of a simple system. In contrast, when considering 
partitioning of many body systems, we are interested in 
scaling properties. Here, the entanglement content of a 
single system grows in a non-trivial way as the system 
size is increased. Thus, we are considering the entropy 
scaling in many-body systems, when super-selection rules 
are applied. 



DERIVATION OF THE MAIN RESULT 

We now proceed to introduce the super-selection ac- 
cessible entanglement entropy. We assume that the local 
operations we may perform in A and B respect conser- 
vation of quantities Na and Nb- The super-selection 
sectors are conveniently represented as the invariant sub- 
spaces of a an operator N = Na <& I + I ® Nb (where / 
is the identity operator). Thus N represents the globally 
conserved quantity (i.e. total number of particles) and 
Na (Nb) are the locally conserved quantities, which can 
be thought of as "number of particles in A (B)" . The 
available entanglement entropy is obtained by averaging 
over the entropy obtained by first restricting the density 
matrix to a sector with fixed numbers of particles N A , 
Nb [9(. More concretely, this is done using the density 
matrices 



1 



Pn.m — 

Pn.rn 



Pn.r. 



nf 



TrTL>n^TL>n^), 



(7) 

(8) 
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where II„ are projectors onto sectors with fixed particle 
number n in A (i.e. all states ip in A, such that Na^P = 
nip) and similarly 11^ projects on sectors with iVg = m 
in B. 

Following Ref.[9|, we define the accessible entropy as 



(9) 



where (/9 n , m )yi = TrBPn.m- In terms of the probabilities 
([8|) we can write the Shannon entropy 



*Sm Pn.m logPn.m 3 



(10) 



which gives the measurement entropy of Na and TVb [23| . 

We now turn to the proof of our main result ([T]). We 
assume that we know that the state of the system \ip) 
is characterized by a fixed total number of particles, 
N\ip) = N\ip). The Schmidt decomposition of the state 
\tp) may be written as \ip) = J2 n K a )A ® \X — n, o)b, 
where a enumerates different states with the same quan- 
tum number, which equals n in A, and N — n in B. Trac- 
ing over B we see that only the reduced density matrices 
{Pn,m)A with to = N — n are nonzero, giving 

S T = - y^Pn,jV-nTr(/>n,A'-w)^ log(p„,jV-n)yl- (H) 
a 

Using the constraint n + m = N, we can also simplify the 
expression ([7]), writing it as 



(ftn,N-n)A 



Pn.N- 



(12) 



Combining this formula with the representation 

PA = J2,iPn,N-n(pr h N-n)A, and noting that 

{pn,N-n)A(pn>,N-n')A = if n ^ n' , we can write 
the accessible entropy © as 



Tr(n n p A n„) log 



Pn,N-n 



(13) 



= -Tip a log Pa + y^P»,jy-n logp„^_„ 

Here 5 ra = - J2n Pn,N~n logp„,Ar-„ is the measurement 
entropy (jlpp for fixed total particle number N = n + to. 
The latter constraint makes particle number measure- 
ments in A and B perfectly correlated, and so particle 
fluctuations in A alone are sufficient to determine S m . 

The relation (|13p may be interpreted as an extension 
of the well known relation from information theory to 
the quantum case: If S(X, Y) is the Shannon entropy 
associated with the joint probability distribution of the 
random variables X, Y, then S(X, Y) = S(X\Y) - S{Y) 
where S(X\Y) is the entropy of X conditioned on know- 
ing Y, and S(Y) is the entropy of Y. 



Next, we estimate by how much S A CS can depart from 
Sa ■ This can be achieved using the relation S^T = 
Sa — Sin derived above, and maximizing S m . It is a basic 
result of information theory that given a variance Ci and 
mean C\, a continuous distribution maximizing the en- 
tropy is a Gaussian distribution (this can be easily proved 
using Lagrange multipliers). The entropy of a Gaussian 
distribution, So = ilog27reC2, where C2 = ((N\)) thus 
supplies an upper bound on measurement entropy of a 
continuous observable. For a discrete variable, such as 
particle number, charge, or spin, this upper bound on 
entropy has to be slightly modified [24[ to be 

S m < AS = 1 log [27re(C 2 + i)] (14) 

establishing the inequalities ([T]). 

It is interesting to note that in Ref. [25| the variance 
of particle number of noninteracting fermions was shown 
to provide a lower bound on the entanglement entropy, 
Sa > (4 log 2)6*2. This inequality can be used to write 



Sa - ST < I log 



2?re 



1 



4 log 2 12 



(15) 



S A CB becomes 



which indicates that the difference Sa 
small in a relative sense as Sa increases. 

It is straightforward to generalize the result ([T]) to the 
case of several conserved quantities, which we refer to 
as "charges" a\ ... a^. For the measurement entropy 
S m (ai...ak) of these quantities, using the subadditivity 
property of entropy, we have 

Sm(ai...a k ) < S m (ai) + ... + 5 m (a fe ). 



Using for each of the terms the inequality (| 14[> . we arrive 
at 

SaST = S m ( ai ...a k ) < [ 2 ^ e ( C M) + M > 

i=l...k 

where 62(0^) is the variance in measurement of the 
charge a% (i = l...k). 



MEASUREMENT ENTROPY FOR STATIC AND 
DYNAMICAL ENTANGLEMENT GENERATION 

Now, having established the connection between ac- 
cessible entanglement entropy and measurement entropy, 
we turn to discuss properties of the latter. In addition to 
being useful as a tool in providing the above bounds on 
S t a s , the quantity S m is directly measurable, and thus is 
of interest in itself. 

In general, of course, the conserved charge does not 
have to be particle number. For example, the conserved 
quantity may be chosen to be energy. In this case, consid- 
ered in Ref. [26[ , it was argued that for some situations 



5 



of interest the measurement entropy of energy S m can 
mimic the behavior of thermodynamic entropy such as 
the second law of thermodynamics, and invariance under 
adiabatic evolution of the system. 

A useful tool for computing S m in many-body sys- 
tems is provided by the generating function, defined as a 
Fourier transform of the probability distribution: 



X(A) 



EA\n 
Pn,N- n e 



Tr (V A %) 



(16) 



Since the generating function p^|) is represented as an 
expectation value, x(A) = {e lXNA ) p , it is often more easy 
to evaluate the quantity (fH))) than the probability distri- 
bution itself. For fermionic systems, in particular, this 
can be done with the help of the determinant representa- 
tion used to analyze counting statistics of current fluctu- 
ations [i^ ■ Once the quantity x(A) is known, the prob- 
abilities p n ,N~n, found from its Fourier transform, can 
be used to evaluate the measurement entropy S m . How- 
ever, if the generating function happens to be Gaussian, 
X(A) = e~° 2X / 2 , the quantity S m can be evaluated di- 
rectly as S m — I log(27reC2) to leading order in C 2 . 

Below we apply this approach to several cases of inter- 
est. We first consider the correction due to S m for for 
entanglement entropy generated in the process of con- 
necting two initially separate systems [22| . The entangle- 
ment generated in this case is logarithmic in observation 
time. In another case, considered in [l2T ] entanglement is 
generated per unit time in a nonequilibrium, but steady 
state of an open system. Here particles are transmitted 
across a QPC in the presence of a bias voltage. In [l2| 
the role of super-selection by a stronger constraint, that 
of energy conservation, was computed, leading to a cor- 
rection to entanglement entropy also scaling linearly with 
measurement time. We find that if only number conser- 
vation is taken into account, the correction to accessible 
entropy is in fact sub-leading. 

We then proceed to consider the accessible entropy of 
critical systems. Such systems are of particular interest in 
the theory of entanglement entropy due to the presence of 
enhanced, non local, correlations 

Here we focus on an interacting system in Id, namely 
a Luttinger liquid, whose entropy can be obtained us- 
ing conformal field theory methods. Somewhat surpris- 
ingly, the accessible entropy of this strongly correlated 
system can also be found since S m may be evaluated us- 
ing bosonization. Finally, we consider a gapless system 
of free fermions in higher dimensions. We show that, in 
an arbitrary space dimension d, the quantity S m can be 
estimated using Widom's conjecture [36| . In all of the 
cases, S m turns out to be a sub-leading contribution to 
the scaling of many-body entanglement entropy. 

(1) First, we consider two Fermi seas coupled through 
a QPC with an externally controled transmission coeffi- 
cient (Fig.[T]a,). Initially the two Fermi seas are discon- 
nected and there is no entanglement between them. The 



QPC is then opened during the time interval < t < At, 
and then closed again. The resulting entanglement en- 
tropy was computed in [22 1. 

The measurement entropy of particle number fluctu- 
ations the Fermi seas, which gives the correction to en- 
tanglement entropy due to charge conservation, can be 
evaluated as follows. We use the generating function 
x(A) defines by (Tl6|) . where the state p is the density 
matrix of the two-lead system at the final time t = At, 
and Na = J x< q ^ (x)^p(x)dx is the charge operator in 
the left lead. The quantity xW m this situation was 
analyzed in Ref.[13| with the help of an expression [3l| 



A/ 



X(A) = exp 



A/logf 
2tt 2 



A* 



D sin ^,(17) 



where D is the QPC transmission during < t < At, and 
t is a short time cutoff of the order of the time it takes 
to switch the QPC from a disconnected to a connected 
state. In this case C2 — -§1 log — . For large At, the main 
contribution comes from small A, and the distribution is 
Gaussian to a good approximation, in formal analogy to 
the Central Limit Theorem. 

We therefore conclude that S m is given by the en- 
tropy of a Gaussian probability distribution, i.e. Sm ~ 
ilog2-7reC2 ~ \ log log Since the full entanglement 
entropy generated in this processes is S ~ \Ci 22], we 
see that the S res = S to leading order in C2. This is our 
first concrete example for our general expectation that 
super-selection leads to a small correction to entangle- 
ment entropy. 

(2) Next, we analyze entropy generated by a voltage- 
biased QPC (see Fig. [!})). To leading order in the time 
of measurement At, we can characterize the systems by 
N = AteV/h independent transmission attempts, with 
the probability of success equal D for each attempt. (The 
effects of QPC opening and closing are sub-leading since 
they are logarithmic in At, as in the previous case.) The 
probability to transmit n particles out of N attempts is 
given by 



Pn 



(N -n)M 



D n {\ - D 



,N-r 



(18) 



The large-iV asymptotic form of the measurement en- 
tropy S m = — ^p„logp„ of overall transmitted charge 
is given by [32j 



Sm = ^log(2neD{l -D)N), 



(19) 



whereas the entanglement entropy, generated in this pro- 
cess is 121 



S = N{D\ogD + (!-£>) log(l - D)). 



(20) 



We see again that S m , Eq. p^l) . is a sub-leading correction 
to the entropy (|2T))) . 
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One may further restrict the allowed operations by not 
considering entanglement between particles that cross 
the barrier with different energies. This restriction, 
which was considered in [l^], substantially reduces the 
entropy, since the measurement entropy of energy turns 
out to be also linear in N. 

(3) Now we shall analyze an interacting fermion sys- 
tems in one dimension, described by Luttinger liquid 
model. Below we shall evaluate the correction to the 
entanglement entropy of a region A of length L due to 
fermion number conservation. In bosonization frame- 
work, the chiral Luttinger liquid is a conformal the- 
ory with central charge c = 1. Therefore its entangle- 
ment entropy scales as |log(L/a), where a is a short- 
distance cutoff 3. However, since the Luttinger system 
is fermionic, this entanglement entropy again includes 
sectors which mix number of particles. Since the overall 
charge of the system is still a good quantum number, the 
relation (fl3|) holds, and so we can estimate the differ- 
ence Sa — S A S by evaluating S m for this situation. The 
Luttinger Hamiltonian can be written as 33] 



H = 



dx: -d x Q{xf + gd x ^(xf 



where Q(x) and d(f>{x) are canonical conjugates and g 
depends on the interaction. Within bosonization ip{x) oc 
E„odd e_iv/ ^ 0(;E)e_m(v/ ^ e(:l;)+fcF;r) - The density is given 
by p = d x Q(x) + k F /n and so N A = e(l)-e(0) + k F L/n. 

The generating function p6[) can now be computed 
from the bosonic theory: 



XLuttW 



= (e iXNA 



iXkpL 1 7r- 



(21) 



giving a Gaussian distribution of charge. Therefore S m 
for this case scales as loglogfcpi 

(4) The last case we shall consider is a problem of free 
fermions in d dimensions. Here the state of the system is 
described by a Fermi sea 



11410), 



ker 



where L is a domain in momentum space, illustrated in 
Fig. [21 defining the set of occupied states. The entangle- 
ment entropy of a region A in real space has been studied 
111 [ItJ in the asymptotic limit where the linear size 
of A is rescaled by a large factor L. 

To compute the correction due to charge conservation 
we again evaluate the measurement entropy S m using 
x(A). The asymptotic form of the generating function 
x(A) at large L can be analyzed as follows. Let Pla be 
a projection on the region A in real space rescaled by a 
factor L (i.e. the set of points Lx where x S A, as il- 
lustrated in Fig|2]). Let Pr be a projection on the set T 
in momentum space, obtained from the Fermi distribu- 
tion (f2"2")l . The generating function x(A) can be written 



Region in real space 



Filled Fermi sea 



LA 



La 



FIG. 2: The region LA, in which particle fluctuations are 
analyzed, is shown along with the region A of the same shape, 
having size of order unity. The scaling factor L is used to 
define scaling of particle number mean and variance, Eq. (|25[l . 
A region F in momentum space defines the Fermi sea, Eq. (|22[) 



as 34, 3 



logx(A) = log(e tXNLA ) = Trlog(l - P r + P r e 
= Trlog(l + P LA P T P LA {e iX - 1)) 



: XPl 



(23) 



where N LA = J LA ■ 

We now estimate x using Widom's conjecture [36j . 
This method has been used in 17J to estimate the scaling 
of entanglement of free fermions in arbitrary dimensions. 
While a rigorous proof for the conjecture is still missing, 
it seems to be perfectly consistent with numerical com- 
3£ 



putations 
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Widom's conjecture [36| is a general- 
ization of the strong Szego theorem to higher dimensions. 
It states that given a function f(z), which is analytic on 
\z\ < 1, with f(z) = 0, the following holds as L — * 00: 



Tr f (P r P LA Pr) = c 1 f(l)L d + c 2 U{f )L d - 1 log L (24) 

f 1 /(*)-*/(! ) 
/o *(1 - t) 



o{L d -"\ogL), U(f) 



■At, 



Here the notation g — o(h) means that g/h — > when 
L — > 00, and the coefficients Ci t 2 are given by ci = 
hfr dxd P> c 2 = jt^rr J dA J dr \n x ■ Up\dS x dS p , 
where n^, n p are unit normals to dA, dT, respectively. 
The formula was proved for d = 1 in [3a ] . 

Plugging f(z) = log(l + z(e lX - 1)) in (EJ), we find 



(22) that U(f) = and so 



logx(A) = i\ Cl L d - yc^- 1 logL. 



(25) 



Thus the charge distribution is Gaussian to leading order, 
and S m scales as log(L d_1 logi) and is again a sub- 
leading correction to the entropy. 

In summary, particle number and charge conservation, 
as well as spin conservation, under some conditions, is 
an essential part of realistic Fermi systems. However, 
the existing discussions of many body entanglement of 



fermions ly, 17 



22| implicitly assume that such conser- 
vation laws have no direct effect on the scaling of many- 
body entanglement. The analysis of the measurement 
entropy of particle number fluctuations and of its rela- 
tion to super-selection rules and entanglement entropy, 



7 



presented above, indicates that this expectation is in fact 
correct. We analyzed several systems of interest, includ- 
ing time dependent scattering problems, the ground state 
of one-dimensional interacting fermions (a Luttinger liq- 
uid), and a free fermion system in higher dimensions. In 
all those cases we found that super-selection rules yield 
sub-leading corrections to the entanglement entropy. We 
conclude that for a generic Fermi system in dimension 
d > 1 the accessible entropy is equal, to leading order in 
system size, to the full entanglement entropy. The same 
is true for critical systems in d = 1 described by con- 
formal field theories with density correlator decaying as 
1/r 2 . 
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